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Abstract. We propose a method of testing source evolu- 
tion theories that is independent of the effects of inhomo- 
geneity, and thus complementary to other studies of evolu- 
tion. It is suitable for large scale sky surveys, and the new 
generation of large telescopes. In an earlier paper it was 
shown that basic cosmological observations — luminosity 
versus redshift, area distance versus redshift and number 
counts versus redshift — cannot separate the effects of 
cosmic inhomogeneity, cosmic evolution and source evo- 
lution. We here investigate multicolour observations, and 
show that by comparing luminosity versus redshift in two 
or more colours, contraints can be placed on source evolu- 
tion even if unknown source evolution is present, providing 
an important test of evolution theories that is complemen- 
tary to present methods. However, number counts in dif- 
ferent colours versus redshift are not useful in separating 
the effects of source evolution and inhomogeneity. 
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1. Introduction 

As measurements of the Cosmos become more extensive 
and more accurate, it becomes increasingly important to 
take account of the evolution of galaxies and clusters of 
galaxies, since these have a direct effect on our cosmo- 
logical measurements. It seems that the Hubble constant 
may soon be well known, and limits on the acceleration 
parameter may improve rapidly. The next step is to esti- 
mate deviations from homogeneity on the large scaled]. It 
is expected that the Sloan digital sky survey (SDSS) will 
detail the galaxy distribution out to z < 0.5. How does 
this data relate to the spacetime metric? 

It was previously shown by Mustapha Hellaby and El- 
lis ( 1997 ) (MHE) that observations of luminosity versus 
redshift and area distance (or luminosity distance) versus 



For example, Wang, Spergel and Turner ( 1998 ) have esti- 
mated that variations of up to 10% between local measures of 
Ho and its large scale value are possible. 



redshift cannot distinguish between the effects of source 
evolution, cosmic inhomogeneity, and cosmic evolution. 
Riess et al Ql998| ) and Perlmutter et al ( |1999| ) have fitted 
supernova observations to an FLRW model with non-zero 



A, but Celerier (2000) has shown they could equally well 
be fitted by an inhomogeneous model with zero AQ. 

Whilst the Copernican Principle may lead many to as- 
sert that homogeneity is valid on a large enough scale, it 
is better not to assume something if it can be measured. 
In addition, the observed luminosity versus redshift re- 
lation does not fit the FLRW predictions without some 
adjustments, such as a cosmological constant that has a 
problematic value, or source evolution. And since source 
evolution is really not known, one should be suspicious of 
the results. Thus it has never been observationally deter- 
mined on which scale homogeneity is valid^|. 

Galaxy evolution is now a very active field, but there 
is a long way to gcQ. There are many functions to be de- 
termined — the rate of formation of each galaxy type, the 
rate of formation of stars of different masses in each galaxy 
type the rate of galaxy mergers^], the effect of galaxy 
mergers and encounters on star formation, how central 
bulges and bars form, how low surface brightness galaxies 
fit into the picture, etc. The early appearance of galaxies 
and the relationship between quasars and galaxies is really 
not known. Many studies of source evolution make an as- 
sumption of homogeneity at some point, by using observa- 



2 Also Maor, Brustein and Steinhardt (2000) argue these 
observations cannot determine the fut ure fat e of the universe. 

3 For example, Barrett & Clarkson (1999) have shown that 



Stephani models with significant inhomogeneity can fit a range 
of key cosmological observations. 
4 see e.g. Ellis, Abraham, Brinchmann and Menanteau 



Impey an d Both un ( p.997|) , Bell and Bower (pOOCh, Bh r 
Rawlings ( |l99S| ), Brunner Szalay and Con nolly ( J2000|) 



fl2000D and referenc es th erein, Kennicut (|1998|) , Ellis ( p,997| ), 

Blundell and 
Bullock 

Dekel and Prim ack ( 200C ) , Bunke r et al ( 200C ) , Cava liere and 
Vittorini ( |200C| ), D ickin son ( |2000| ), Fasano et al (|200C| ), Genzel 
Lutz and Tacconi (199£ ), K odama and Bower ( 2000|), Ponman 
Cannon and Navarro ( |1999| ), Watts and Taylor ( ]200c| ). 

5 It has also been proposed that stellar mergers may be 
significant in dense clusters. 
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tional relations derived in a FLRW model. While we don't 
dispute their usefulness, the danger is that the resulting 
source evolution relations may later be used to demon- 
strate homogeneity. Thus methods of determining source 
evolution that don't make assumptions about the universe 
model should be emphasised as more reliable. This point 



was also made by Goodman (1995). Observations of super- 
novae (Riess et al 1998, Perlmutter et al 1999) may help to 



separate inhomogeneity from source evolution, out to dis- 
tances where their light curves can be measured, though it 
is still uncertain whet her they are free of evolution effects 
(e.g. Riess et al |l999[) . 

As a contribution towards this requirement, we inves- 
tigate whether multicolour observations can help to sepa- 
rate inhomogeneity from source evolution. The measure- 
ments in different colours may be a selection of spectral 
line intensities or luminosties in U, B & V filters, say. In 
any case, one needs to know the absolute luminosity in 
each colour at each z value, so there is a new source evo- 
lution function for each colour. Thus it seems that multi- 
colour observations may not improve the situation. How- 
ever, a key point is that the luminosity distance must be 
the same in all wavelengths for each given source. So if the 
luminosity-redshift plot in two colours is not the same, this 
is evidence that the two colours have different evolution 
functions. This approach is used below, and it is found 
that the basic uncertainty in absolute values remains, but 
the luminosity-redshift relations in different colours must 
be related to each other, thus providing constraints on 
evolution theories. 

The choice of cosmological model is not central to the 
ideas presented here, but to keep the equations simple and 
focus on the basic concepts, we choose the simplest inho- 
mogeneous cosmological solution of Einstein's equations. 
This is the Lemaitre-Tolman (LT) model^j — a spherically 
symmetric dust cosmology that is inhomogeneous in the 
radial direction. If one accepts that the large scale universe 
is a collection of galaxies with negligible bulk rotation and 
interacting only through gravity, then the dust equation 
of state is valid, and the cosmic time evolution is pretty 
much determined along each worldline. For reasons of sim- 
plicity we assume the observer is at the centre; thus we 
are assuming large scale isotropy about the milky way or 
some nearby point. Isotropy is relatively easy to test and 
is not a bad assumption on large angular scales. It is also 
a natural simplification to make, since we are at the cen- 
tre of the null cone we observe. A more general approach 
would cloud the issue at hand, though it will be needed 
in the long run. More importantly, the degree of deviation 
from isotropy is directly seen and is not mixed up with 
source evolution or cosmic evolution. 



Though it is sometimes calle d the Lemaitre-Tolman-Bondi 
(LTB) model, Krasinski's (1997) terminology is adopted here. 



2. The Cosmological Model 

We summarise the results of MHE, as background to the 
present considerations, and refer the reader to the papers 
cited there. 



The Lemaitre-Tolman (LT) metric (Lemaitre 1933 
Tolman |1934 Bondi |1947|) is 



ds 2 = -dt 2 



{R'f 



■ dr 2 + R 2 dfl 2 



(1) 



1 + 2E 

where R = R(t,r), E = E(r), R' = dR/dr, and dVL 2 = 
d6 2 + sin 2 9 dtp 2 . The function R is the areal radius, which 
obeys 

2M 
~R 

and the arbitrary function E determines both the local 
geometry and the local type of time evolution. The ar- 
bitrary function M = M(r) plays the role of the total 
gravitational mass interior to shells of coordinate radius 
r. This has solutions in terms of the parameter r\ 
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and ts = £b(?") is a third arbitrary function. It gives the 
time of the big bang, R = 0, locally. Although the initial 
singularity is spacelike everywhere^ the spacetime emerges 
from the initial singularity over a finite or possibly infinite 
time. The density is 



8np - 



2M' 
R 2 R' 



(5) 



Putting the observer at the centre of symmetry, r — 0, the 
observer's past null cone is the solution of 



dt 



-R' 



VTT2E 



■ dr 



(6) 



that passes through r = at t = to — i.e. at the present. 
MHE wrote this particular solution as t — t(r), and a hat 
will henceforth indicate quantities evaluated on this null 
cone. Since we really only need this one path, MHE used 
the freedom in the radial coordinate r to specify 



R' = l 



(7) 



on t — t(r) only. This greatly simplifies the equations to 
be solved. The path of the light cone is then 



f (r) = t - r 



(8) 



except where there's a shell-focussing singularity 
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MHE then went on to show how the 3 arbitrary functions 
of the LT model, E(r), M(r) and t B {r), could be deter- 
mined from observations, if one knew the source evolution. 
The observations needed were the source number counts 
against redshift z, and the diameter distance R (or the re- 
lated luminosity distance) against z. The evolution func- 
tions needed were the mass per source, and the absolute 
diameter or luminosity, both against z. Given only these 
observations, it was shown that both source evolution and 
inhomogeneity cause deviations from the expected FLRW 
observational relations^, and there is no way to distinguish 
the two. 

We now generalise the MHE approach to observations 
in two or more different colours. 



3. Observing Model 

3.1. Assumptions 

We assume the following: 

— The universe is isotropic about the earth to good ap- 
proximation on the large scale. 

— Small scale inhomogeneity has been eliminated by av- 
eraging cosmic observables over the whole sky, so that 
they are functions of z only. 

— There are J types of visible sources, 1 < j < J , 
e.g. spirals, ellipticals, Abell clusters, field galaxies, 
etc. The various types can be reliably distinguished. 
The types have different intrinsic properties and may 
evolve differently — i.e. have different source evolution 
functions. Mergers could correspond to the removal of 
certain types in favour of new types. 

— There are I different spectral frequencies or colour fil- 
ters being used for observing, 1 < i < I , e.g. U, B & 
V filters, or preferably a set of spectral line intensities. 

— All wavelengths experience the same delay and the 
same cosmological redshift due to the geometry and 
evolution of the universe^. 

3.2. Notation 

Whether the colour measurements are spectral line inten- 
sities or colour filter apparent magnitudes, we will call 
them "colours". 

3.2.1. Observables 

The following "observables" are smooth functions fitted to 
data from direct observations, that have been corrected for 



8 Examples of strong deviations from the standard FLRW 
observational relations in models with realistic amounts of in- 
homog eneity were given in Mustapha Bassett Hellaby & Ellis 
( [l998| ) (MBHE). 

However redshifts of individual objects due to peculiar 
motions will differ, as noted below. 



selection effects, absorbtion and other effects mentioned 
below. 

— z — measured redshift of each object. 

— nj{z) — observed number density of source type j in 
redshift space, per steradian per unit redshift interval, 
1 < j < J ■ 

— (ji(z) = measured apparent luminosity of source type 
j in colour i at redshift z, 1 < i < J . 

— 5ji(z) = measured angular diameter of source type j 
in colour i at redshift z. 



3.2.2. Theoretical Quantities 

The following are functions that must be supplied by a 
theory of source evolution. 

— Lj(z) = absolute bolmetric luminosity of source type 
j, at the epoch corresponding to redshift z. 



Li 



absolute luminosity of source type j, in colour 



i, at z. 

— Dji(z) = proper diameter of source type j, in colour i, 
at z. 

— mj(z) — total mass associated with (gravitationally 
bound to) source type j. 

— v(z) — total proper density in redshift space of matter 
not associated with a luminous source. 

3.3. Complications 

There are of course many sources of error in making and 
reducing the observations, which need careful attention. 
However we will not go into them here, except to note 
the most important ones. It is assumed the observational 
functions have already been corrected for these effects. It is 
evident that measurements of spectral line intensities have 
fewer problems than luminosities through colour filters. If 
a program of observations were being planned, line inten- 
sities are a natural extension to redshift measurements, 
and would be the natural choice. However measurements 
of colour magnitudes are considerably easier to make, and 
more such data already exists. 

— The redshift we are interested in is the cosmological 
one, but the actual redshift of a source is the combi- 
nation of cosmological and peculiar velocity contribu- 
tions, causing considerable scatter in the observed z 
values from a given distance. Since the functions of z 
we use here are all-sky averages, this is not a problem. 
However, in generalising to models with a measure of 
anisotropy in addition to radial inhomogeneity, a net 
'peculiar velocity' on the averaging scale being consid- 
ered, should rather be viewed as an inhomogeneity in 
the cosmological expansion rate^ 



10 However, as was shown in MBHE, inhomogeneities near 
the maximum in R(z) can create loops in the R(z) graph, so 



4 



Charles Hellaby: Multicolour Observations, Inhomogeneity & Evolution 



The observed colours will be redshifted. If spectral line 
intensities are used, then the emitted frequencies will 
be known, but if colour filters are used, different parts 
of the true spectrum are selected at each z value, so 
correcting for this could be quite tricky, to say the 
leastQ- 

Sources will suffer absorbtion and reddening. This 
could be problematic as intergalactic absorbtion is not 
all that well known. Ideally wavelengths should be cho- 
sen to minimise absorbtion. 

Selection effects ar e cruc ial to real obser vations (El- 
lis, Per ry fc Sievers 1984 , Teerikorpi 1997 , Totani and 
Yoshii |2000| ). Only a fraction / of sources are actu- 
ally detected, and / depends on a source's apparent 
size and surface brightness. Thus the observed number 
densities rij are related to the true number fij densities 

by 



(9) 



For quantities calculated from the rij, this is a serious 
problem. But for those that depend on comparisons of 
measured Iji values — exactly what we need below — 
there is no real problem. 

Source evolution is naturally a function of time at the 
source, r, but here it is written as a function of z. Al- 
though the age of the universe at each z is an output 
function of the MHE algorithm, one needs the source 
evolution functions as initial inputs. An iterative cor- 
rection process could be developed, once our knowl- 
edge of source evolution is adequate to the task. 
It could be difficult to identify the same source type at 
different stages of evolution, especially at high z values, 
and especially if galaxy mergers are very frequent. For 
this reason it may be best to initially consider a single 
source type, and test for the bulk evolution properties 
of the cosmic luminous matter content. 



4. Comparison of Multicolour Observations 

Multicolour observations cannot be directly used as input 
to the LT model (or any other metric) , which only requires 
R(z) and p(z) (plus the coordinate condition R' — 1) to 
fully determine the metric. However for any given source, 
the ratio of apparent to absolute luminosity must be the 
same in all wavelengths. In other words, the luminosity 
distances obtained in each of the colours or spectral lines 
must be the same. 



that objects at different distances have the same z. This is 
very like the "finger of god" effect. Thus, if this were common, 
sources selected by z may not all be at the same true distance, 
and there would be a blurring of source properties near the 
maximum, in addition to any observational uncertainties. 

However, a well developed source evolution theory could 
make predicions for the relevant part of the spectrum according 
to z. 



4-1. Distances 

The luminosity distance is 



d L = 



L 

~4n£ 



(10) 



Ideally, all estimates at any given z, for all source types 
and all colours, should agree, if we knew the Lji(z) com- 
pletely: 



(11) 



This is a set of / x J equations for the I X J + 1 unknown 
functions Lji(z) and (1l(z). Thus with sufficient observa- 
tions we could determine the Lji(z) relative to each other, 
but none of them absolutely. 

To get the function R(z) we take the average^ 



R(z) (1+z) = d L , d L = 



l-i , 



(12) 



In reality, observations would be of different reliability, 
and a weighted average would be used. 

We have similar relations for the diameter distance 



D 



all of which should agree, if we knew the Dji(z) 
Dji(z) 



Sji(z) 



d D (z) 



(13) 



(14) 



Again the average should give a good estimate of the LT 
areal radius 



R(z) = d D 



EJ f l v^J Dji_ 



(15) 



and no absolute diameters can be determined. 

If both types of distance measurements are available, 
then we have a combined requirement: 



1 ^ =R\z)=(^- 



(1 + z) 2 Anljiiz) 



Vi.j (16) 



leaving only one free function out of the set {Lji, Dji}. 

From here on we will focus on the luminosity distance, 
as the one that can be used to the greatest depth. 



4-2. Densities 

If the number of sources of type j observed between z 
z + dz in solid angle dfl is 



dfl dz 



(17) 



Equation (31) of MHE is incorrect. 
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then the net mass is 



5.1. Colour-Weighted Number Counts 



rij rrij dil dz 



(18) 



so the combined mass of all source types seen over the 
whole sky is 



4n± 

j=i 



Uj rrij dz 



(19) 



and if we allow for a distribution of non- visible matter that 
isn't associated with any source, we get the total mass as 



Number counts of field galaxies with different colours have 
been used to try and demonstrate source evolution. If 
such a method is undertaken, it is because the iji are not 
known, and also because spectral line data is not available, 
only colour filter luminosities. We here try to model this 
process in an admittedly rough and ready manner — the 
number counts obtained in different colours are modelled 
by weighting the total number count by a function of all 
the colour luminosities of that source type. 
Let 



4-7T I rij rrij +v J dz 

Now the local proper density in the LT model is 
P = P{t,r) 

and its value on the null cone is 
P = P(i{r),r) ■ 

The total mass between r and r + dr is 

— . ^R 2 R> 
pd 6 V — p - dr 



VT+2E 



and hence 

„ R 2 R> 
9 VI + 2E 



dz 



\ - 1 uz 



(20) 

(21) 
(22) 

(23) 
(24) 



4-3. Fitting the Data with an LT Model 

Equations @ and @ or @ or (JX6j) show how the LT 
functions relate to the observables and the source evolu- 
tion functions. Once the last two are known, the MHE 
procedure shows that, if we are given 



nj(z) mj{z) + u{z) and d(z) 



(25) 



then the 3 arbitrary functions that characterise the best fit 
LT model can be determined. Thus if we knew the source 
evolution, we could determine the LT model that repro- 
duces the observations. In particular, we could determine 
how close the universe is to large scale homogeneity on 
the smoothing scale used. 

5. Number Counts 

The number counts only provide a single function of z for 
each source type, and so do not help us to distinguish in- 
homogeneity from source evolution. The advantage of the 
luminosity measurements is that each source is measured 
more than once in different frequencies. 



— Wji = weighting functions for each source type. These 
may in general depend on all / of the colour luminosi- 
ties, e.g. w 2 3 = ^23(^21,^22,^23), 

= number of source type j, counted in colour band 



J* 



Thus at each z value: 



n ji — w ji(Lji) Tlj 



(26) 



Cj = colour-weighted number count of all source types 
in colour band i — i.e. the number density in redshift 
space Ci{z) of all source types with colour i is 



J J 
4=1 i=i 



Wji(Lji) rij 



(27) 



The simplest weighting function is a step function with a 
cut-off luminosity below which the source is not detected, 



1 L 
L 



(28) 



This means that each of the equals the total number 
count out to some z value, and after that each source type 
disappears from the count fairly rapidly around some z 
value. The z values for each type and each colour are prob- 
ably not very different. Again any variation could equally 
well be due to inhomogeneity in the population of source 
types as to source evolution. Add to this the difficulty of 
correcting for the redshifting of frequencies detected in 
each colour band, and it is clear that this data does not 
assist in distinguishing source evolution from inhomogene- 
ity. Indeed, even if there were no large scale inhomogene- 
ity, it is still difficult to get clear information about source 
evolution from this data|^|. 

Now suppose the weighting function were proportional 
to the relative luminosity of that source type in that colour 



Wji 



it** 



(29) 



For example, morphological studies suggest that the excess 
of faint blue galaxies is due to those with peculiar morphology. 
But which of the more modern nearby galaxy types they evolve 
into, if any, or whether this is actually a manisfestation of 
cosmic inhomogeneity, is not known. 
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This is of course not realistic, but by presenting an exces- 
sively favourable case we highlight how unpromising the 
more realistic scenarios are. For these weighting functions 
we would have 



v 3 l 



L 



1 



7 I> 

i=l 



(30) 



and so the number density in redshift space Ci{z) of all 
source types with colour i would be 



j 

n = E 



-=d-n. 



(31) 



where Jj and Lj are the mean luminosities. Assuming the 
rij are known, we have / linear relations involving I x J 
unknown coefficients. This certainly doesn't help us solve 
for the evolution functions, though it does place mild con- 
straints on them. The ideal situation is if there is only 
one source type, in which case we have I equations for I 
unknowns. 

Any realistic weighting function would smear together 
these equations, making it virtually impossible to solve for 
the the coefficients with any certainty. 

These results, even for the most ideal case, are condi- 
tional on the number counts in each colour being tabulated 
against redshift z. But since they are actually summed over 
a large range of z, then the constraints are lost, as this ef- 
fectively smooths out the evolution functions over time. 

6. Examples 

6.1. 1 Type of Source & 3 Colours 

This is the case where we treat all galaxies as more or less 
the same, with roughly the same luminosity and colour 
evolution. 

We have J = j = 1, I = 3, i = 1, 2, 3, so all subscripts 
are colour subscripts. The constraints are 



Li 



= d L {z) 



Anil 47T^ 2 47r4 

where all the Ls and £s are functions of z, and 



dL = l 



4tt4 



L-2 



4tt4 



(32) 



(33) 



This gives effectively 2 constraints on 3 evolution functions 
— e.g. L 1 /L 2 and L 2 /L 3 . 

6. 2. 2 Types of Source & 1 Colour 

In this case we can distinguish say 2 types of source, with 
different evolutions, but we only measure total luminosi- 
ties. 



We have J = 2, j = 1, 2, I = 1 = i, so all indices are 
source type indices. The constraints are 



Antx 4tt£ 2 



= d 2 L (z) 



(34) 



This clearly constrains the relative evolution of the source 
types. 

6.3. Several Types of Source and Several Colours 

If we need to distinguish 2 (or more) source types, and we 
make multicolour observations, then we get many more 
constraints, and we improve our chances of demonstrat- 
ing source evolution. If however the evolution of the source 
types is not very different, the extra effort would be un- 
productive. Indeed, the more types we distinguish and the 
more colours we use, the smaller the sample size in each 
redshift interval. However, with the large scale sky surveys 
now operating or being developed, this may not be much 
of a limitation. Nevertheless, identifying different source 
types requires resolved images or detailed spectroscopy 
for morphological classification, and involves considerable 
work. 

7. Conclusions 

— The main result is that measurements of apparent lu- 
minosities in various colours {(ji(z)} — ideally a set 
of spectral line intensities — put strong constraints 
on the colour evolution functions {Lji(z)}. They de- 
termine all of the relative colour evolution functions 
Lji(z)/L(z), but not L(z). The same applies to the 
apparent diameters. This allows source evolution the- 
ories to be tested against the observational data with- 
out any assumptions about homogeneity, or the cosmic 
equation of state. 

— The essential point is, although luminosity distances 
and diameter distances are model dependent, affected 
by both inhomogeneity and equation of state, their 
ratios in different colours are not. 

— The converse is not true — multicolour observations 
do not directly help to pin down the degree of inhomo- 
geneity. Once the absolute source evolution functions 
are reliably known from a well-confirmed galaxy evo- 
lution theory, a fit of the observational data would be 
possible. Otherwise, we would still need to know the 
true distances, by an independent method, such as su- 
pernova light curves, if they prove reliable, or gravita- 
tional microlcnsing. 

— Colour band number counts C{(z) by themselves (with- 
out luminosity measurements) do not put much con- 
straint on candidate evolution theories, even if they are 
known versus redshift z. But, if they are summed over 
a range of z values as is usual they tell us nothing. 

— This method gives no direct information about the 
evolution quantities rrij(z) and v(z). However they 
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are more amenable to determination through measure- 
ments of orbital velocities in galaxies, galaxy interac- 
tions, gravitational lensing surveys, etc. They will also 
be part of the evolution theories that are being tested 
by this method. 

— The background model (LT) is not central to the ideas 
presented here. The idea of comparing observational 
relations in different colours or frequencies to obtain 
clear evidence of source evolution will obviously apply 
in any model. 

— As already emphasised, it is important not to test 
source evolution theories by fitting to a homogeneous 
FLRW model, because then the possibility of detect- 
ing large scale inhomogeneity — or of demonstrating 
homogeneity — is removed, and the theories are wrong 
if there really is inhomogeneity. 

— The approach suggested here offers the possibility of 
testing evolution theories independently of whatever in- 
homogeneity may be present, which is a distinct advan- 
tage over other methods. However, the determination 
of the inhomogeneity remains dependent on knowing 
the source evolution functions quite well. 

— It is worth emphasising that obtaining non-constant 
functions Lji(z)/L(z) from multicolour observations 
does not itself prove the variation is due to time evolu- 
tion — it could equally well be due to spatial variation. 
What it does do is provide data for testing theories of 
source evolution that is not contaminated by the grav- 
itational effects of inhomogeneity in the intervening 
space. 

— The ideal observations envisaged — a redshift value 
and at least two spectral line intensities on a large 
number of galaxies — may be possible, with large 
scale observing programs along the lines of the SDSS, 
or with the new generation of telescopes: Keck, Sub- 
aru, Gemini, VLT. The SDSS0 is compiling an ex- 
tensive and uniform database of galaxy data, for use 
in mapping the galaxy distribution, studying galaxy 
and quasar evolution and luminosity functions, im- 
proving the values of fundamental cosmological pa- 
rameters, constraining dark matter distribution, and 
gravitational lensing studies. When completed, it will 
have photometrically recorded a quarter of the sky in 5 
colour bands between 3000 & 10 000 angstroms. From 
this, 900 000 galaxies with mean z » 0.1 plus 100 000 
luminous red galaxies with mean z « 0.5, and 100 000 
quasars, will have been selected for spectroscopic imag- 
ing. 

— We advocate an observational program, possibly us- 
ing SDSS data, if suitable, aimed at extracting the 
source evolution functions, and ultimately establishing 
the degree of large scale inhomogeneity, through the 
approach suggested here, as a very worthwhile comple- 
ment to other methods being pursued. While there are 



plenty of practical problems involved, they are all fa- 
miliar or being faced in projects currently under study. 
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